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Abstract 

We discuss a neutrino mass matrix M u originally found by Babu, Ma, and Valle 
(BMV) and show that this mass matrix can be characterized by a simple algebraic 
relation. From this relation it follows that atmospheric neutrino mixing is exactly 
maximal while at the same time an arbitrary mixing angle #13 of the lepton mixing 
matrix U is allowed and — in the usual phase convention — CP violation in mixing is 
maximal; moreover, neither the neutrino mass spectrum nor the solar mixing angle 
are restricted. We put forward a seesaw extension of the Standard Model, with three 
right-handed neutrinos and three Higgs doublets, where the family lepton numbers 
are softly broken by the Majorana mass terms of the right-handed neutrino singlets 
and the BMV mass matrix results from a non-standard CP symmetry. 
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1 Introduction 



The atmospheric neutrino problem, with mixing angle #23, requires sin 2 2# 2 3 > 0.92 at 90% 
CL, with a best fit value sin 2 2# 23 = 1, i.e. maximal mixing pQ. There are many models 
and textures in the literature which attempt to explain large — not necessarily maximal — 
atmospheric neutrino mixing — for reviews see Ref. But, the closer the experimental 
lower bound on sin 2 2023 comes to 1, the more urgent it becomes to find a rationale for 
maximal atmospheric mixing. Unfortunately this is not an easy task. Maximal mixing 
means {U^l = \U T s\, where U is the lepton mixing matrix, and this in general requires a 
/i-r interchange symmetry, which on the other hand must be broken since 7^ m T . For 
a recent discussion of this point see Ref. [3] . 

Two models for maximal atmospheric mixing have been suggested by us, one of them 
jl] based on lepton-number symmetries softly broken at the seesaw scale, the other one 
jH] based on a discrete symmetry spontaneously broken at the same scale. Both models 
yield an effective mass matrix for the light left-handed neutrinos at the seesaw scale 



M v = y z w , (1) 




where x, y, z, and w are in general complex. The matrix 

H = M* V M V (2) 

then has an eigenvector (0, 1, — 1) T , and therefore the models predict U e3 = besides 
\U/ia\ = I 1 ; they will have to be discarded if \Ues\ is experimentally found to be non-zero. 

A different approach has been suggested by Babu, Ma, and Valle (BMV) |0j. Starting 
from a degenerate neutrino mass matrix at the seesaw scale, and using the renormalization 
group in the context of softly broken supersymmetry with a general slepton mass matrix, 
they have obtained at the weak scale 



M v = r s b , (3) 




where r and s are in general complex while a and b remain real. BMV have found that, 
under some approximations, this M. v yields maximal atmospheric mixing and, further- 
more, an imaginary U e z ("maximal CP violation") in the standard phase convention (to 
be specified shortly). 

It is important to note that the mass matrix in Eq. (JSJ) is not a generalization of the 
one in Eq. (pj. If r and s in Eq. (JHJ) are real, then that mass matrix coincides with the 
one in Eq. with x, y, z, and w real. On the other hand, the mass matrix of Eq. (jHJ) 
in general yields a non-zero U e $ as soon as r and s are complex, while $1} always yields 
U e3 = 0, even when x, y, z, and w are complex. 

In Ref. [7j, Harrison and Scott (HS) suggested that the lepton mixing matrix may 
satisfy 

\Urf\ = \U Tj \ for J = 1,2,3, (4) 
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which is still a particular case of maximal atmospheric mixing, but has the advantage 
of being more general than the extra condition U e3 = 0. HS showed that, if U satisfies 
Eq. (j3J) and the neutrinos are Dirac particles, then U may be parametrized as 

/ Mi U 2 U 3 \ 

U = \ w 1 w 2 w 3 , (5) 

V wi w* w ; ) 

where the Uj (j = 1, 2, 3) are real and non-negative, while the Wj are complex and satisfy 
the orthogonality conditions 

2 Re {wjW* k ) = 5jk — UjUk ■ (6) 

HS also introduced the concept of "/x-r reflection," which they defined as "the combined 
operation of fi—r flavor exchange [...] and CP transformation on the leptonic sector" and 
which is embodied in the mixing matrix of Eq. (J3J). 

It is the purpose of this paper to, firstly, prove that the mass matrix of BMV al- 
ways yields Eq. (J3J) and, as a consequence, exact maximal atmospheric neutrino mixing 
and maximal CP violation. We shall also show that the BMV mass matrix leads to a 
mixing matrix of the form in Eq. (jSJ), even while the neutrinos are Majorana particles. 
Secondly, we shall put forward a model, based on softly broken lepton numbers and on 
the non-standard CP symmetry called "/i-r reflection" by HS, which obtains the BMV 
mass matrix at the seesaw scale — without the need for the renormalization group, for 
super symmetry, or for an extended fermion spectrum like in the original BMV model. In 
this way we conclude that maximal atmospheric mixing is compatible with a non-zero U e3 
and can be obtained in an extension of the Standard Model. 



2 The mass matrix 

Let M. v be a symmetric complex 3x3 matrix, the Majorana mass matrix of the light 
neutrinos, defined by 

£ mass = - v T L C~ l M v v L + H.c. (7) 

[C is the Dirac-Pauli charge conjugation matrix), in the basis where the charged- lepton 
mass matrix is diagonal. The lepton mixing matrix U is the diagonalizing matrix of Ai u , 
defined by 

U T M. U U = rh = diag (mi, m^, rn 3 ) , (8) 
where the masses rrij are real and non-negative. 

Lemma: Suppose U and U' satisfy Eq. (JHJ) and the masses are non-degenerate. Then 
there is a diagonal unitary matrix X such that U' = UX. Furthermore, Xjj is an arbitrary 
phase factor if rrij =0, while Xjj = ±1 for rrij ^ 0. 

Proof: Since both U and U' fulfill Eq. (JHJ), M v = U*mW = U'*mU'\ or 

W*m = mW with W = U' ] U . (9) 
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This equation, together with the non-degeneracy of the masses, forces W to be diagonal, 
i.e. W* = X. It is moreover clear that Wjj is real when rrij is non-zero, Q.E.D. 
The matrix M. v of Eq. (j3J) is characterized by 



(10) 




SM U S = Ml with S = 

Let us write U = (01,02,03) with column vectors Cj. Equation (JBJ means that 

M u Cj = rrijC* . 

Starting from this equation and using Eq. we see that 



My ( SC* 



m. 



(sc*y. 



in 



(12) 



We thus have a second diagonalizing matrix U' = SU*. Using the lemma above we find 
that, if the masses are non-degenerate, 



SU* = UX . 



(13) 



Consequently, Eq. (J3J) holds. 

An alternative proof of Eq. (0J) starts from the observation that the matrix H corre- 
sponding to the M. u of Eq. © has 



H, 



MA' 



H TT and H efl = H* T . 



(14) 
(15) 



As a consequence, 

(H n )„ = {H n ) TT 

for any positive integer n. Using H = Urh 2 W , it follows from Eq. (J15)) for two distinct 
values of n that either there are degenerate neutrinos or Eq. (@J) holds. 
A popular representation of U [H] is given by 



f/ = diag(V Ql , e ta \ e ta *) U 23 U 13 U 12 dmg(l, e t(3 \ j 



(16) 



with 



U. 



23 



u 



13 



u- 



12 



t 1 \ 

C23 S23 , 

V -s 23 c 23 / 

C13 s 13 e- i<5 
1 
-s 13 e i(5 c 13 

/ C12 S12 
-S12 C12 

V 01 



(17) 



(19) 
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The phases atj (j = 1,2,3) are unphysical (unobservable); 5 is the Dirac phase and /3i j2 
are the Majorana phases. Computing the product of matrices one gets 



U23U13U12 



s 13 e 
S23C13 

V s 12 s 23 ~~ c 12 c 23 s 13 el<5 ~~ C 12-S23 ~~ •Sl2C23 s 13 ei ' 5 C23C13 



C12C13 
-S12C23 - Ci 2 S23Sl3e l 



S12C13 
C12C23 - Sl2S23-Sl3e 1 ' 



(20) 



Equation (j^J) applies to this product. From that equation with j = 3 one obtains 

''23 *23 ■ — • (21) 



1 

V2' 



Now we inspect Eq. (j3J) with j = 1, 2. We know experimentally that C12S12 7^ 0, since 
solar neutrinos oscillate [9 j. It follows that 



s 13 cos 5 = 



'.22) 



i.e. either U e3 

iS 



or CP violation is maximal. Conversely, if we require maximal CP 
violation (e lt> = ±z) and maximal atmospheric neutrino mixing ( | C/^*3 1 = 1^3 1) with the 
parameterization of Eq. (J2UJ) . it is easy to see that Eq. (jlj) follows [7j. 

We stress that the mass matrix of Eq. (jSJ) restricts neither the neutrino mass spectrum 
nor the solar neutrino mixing angle Q\2- In the general case cos 5 = also #13 remains free. 
Note that, with Eq. (|21)l. the parameter measured in atmospheric neutrino oscillations is 
sin 2 2# atm = 4 \U^\ 2 (1 - \U^\ 2 ) = 1 - 4. 

Now we want to discuss the relation between the BMV mass matrix and the param- 
eterization of the mixing matrix in Eq. (jSj). We stick to the — experimentally justified — 
assumption that the neutrinos are non-degenerate, and we employ again Eq. (fT3*j) . If 
rrij 0, then we know that the Xjj are either +1 or — 1. If Xjj = +1 then we see from 
Eq. (0 that 



with real Uj. If Xjj 



— 1 then 



W ; 



lUj 
W ; 



(23) 



(24) 



It remains to consider the possibility m» 



= 0. In that case we have Sc* = CjXjj with 
an arbitrary phase factor Xjj. Since the massless case allows rephasing of the Majorana 
neutrino field, one can absorb a factor (Xjj)~ 1 / 2 into that field. It is easy to see that, then, 
Cj assumes the form in Eq. (J23|) . In the case of Eq. (|24|) . we may multiply the physical 
neutrino field by a factor i, thereby passing from Eq. (|24j) to Eq. ()23j) . but also changing 
the sign in front of rrij in Eq. We may also, if needed, multiply the neutrino fields by 
factors —1 so that all three Uj become non-negative. We thus obtain the following result: 
if the mass matrix is of the BMV type (and since neutrinos are non-degenerate) then the 
lepton mixing matrix U is of the form in Eq. (j3J), but the Majorana phase factors r\j may 
be either 1 or i; or, in other words, with a U of the form in Eq. (j3J), the BMV mass matrix 
is diagonalized as U T M. V U = diag^mi, ^|m 2 , i] 2 m 3 ). 
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Let us finally consider the conditions under which U e3 will be zero with a mass matrix 
of the BMV type. Suppose r 2 s* in Eq. (jSJ) is real. Then, 

(a \r\ \r\ \ 
\r\ ±\s\ b \Y with Y = diag (l , e iargr , e" iargr ) . (25) 
|r| 6 ± |s| / 

This means that M. v is essentially identical to the matrix in Eq. (£[]), and we conclude 
that if r 2 s* is real then U e3 = 0. Conversely, let us now suppose that the neutrino masses 
are non-degenerate. Then we know, from Eq. ([TH|) . that Sc 3 = ±c 3 . (This relation, with 
the plus sign, also holds for m 3 = 0, as we have argued in the previous paragraph.) If 
U e3 = this means c 3 = (0 , w 3 , ±Wg) T . Now, from Eq. (fTTj) . M. u c 3 = 777.303. This gives 

rw 3 ±(rw 3 )* = 0, (26) 
sw 3 ± 6^3 = m 3 ?x>3 . (27) 

Equation (|27]) implies that sw^ is real. Equation ()26)1 implies that r 2 w 3 is real. As w 3 7^ 0, 
we conclude that provided the neutrinos are non-degenerate, U e3 = implies a real r 2 s*. 
We have thus shown that, with the mass matrix of BMV, U e3 being zero is equivalent to 
r 2 s* being real. 



3 A model 

We now want to produce a model that leads to the mass matrix of Eq. Q. In doing 
this we find inspiration in our model of maximal atmospheric neutrino mixing of Ref. j3j. 
Thus, we supplement the Standard electroweak Model with three right-handed neutrinos 
and two extra Higgs doublets. We denote the three lepton families e, /1, and r by the 
general index a; thus, we have three left-handed lepton doublets 



D r 




a: 



e,/i,T. 



:28) 



together with three right-handed charged-lepton singlets and three right-handed neu- 
trino singlets v a R. In the scalar sector we employ three Higgs doublets 



J = 1,2,3 



(29) 



These Higgs doublets acquire vacuum expectation values (VEVs 
and v = 











V 3 



72, 



1 2 . 
V2\ + 



|2 

V 3 \ ~- 



246 GeV represents the Fermi scale. 
We introduce the three U(l) lepton-number symmetries L a . These symmetries are 
meant to be broken only softly at the high seesaw scale. We also introduce a Z 2 sym- 
metry under which fin, tr, 02, and (fi 3 change sign. This symmetry Z 2 is broken only 
spontaneously by the VEVs v 2 and v 3 . Because of the lepton-number symmetries and of 
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the 7L 2 symmetry, the Yukawa Lagrangian of the leptons is (see also Refs. |UE]) 



£y = -— ( Vl, ~(pi ) Y. fcPaRDa ~ ( <Pl ' fT ) C R D e 

3 

£ ( ^7> V?* ) 9j a a R D a , (30) 

where the three f a and the four gj a are complex numbers (m e is real without loss of 
generality and represents the electron mass). Notice that, through the first line of Eq. (|30p. 
the smallness of the neutrino masses may be correlated with the smallness of the electron 
mass. The 1>2 above is analogous to the auxiliary Z2 of Refs. jU E] . 
The right-handed neutrinos have Majorana mass terms given by 

= ~ (f R C- x M* R v R - v R M R Cv T R ) , (31) 

where M R is a 3 x 3 symmetric matrix in flavor space. Now, M R is not diagonal since the 
terms in Eq. (JBTj) have dimension three and we allow the lepton-number symmetries L a 
to be broken softly [3] . Indeed, M R is the sole source of lepton mixing in this framework. 



According to the seesaw formula [TO] , when the eigenvalues of ^M R M R are all of order 
m R ^> v one has 

M v = -MlM R l M D , (32) 

where Md = diag (f e , f^, f T ). It has been shown in Ref. ^1] that this framework, in 
which the tree- level Yukawa couplings are diagonal but M R is not, leads at the one- loop 
level to a renormalized theory with flavor- changing neutral Yukawa interactions, in which 
flavor-changing processes like ^ — > e ± 7 or Z° — > e ± / u =F are suppressed by inverse powers 
of m R while processes like ^ —>■ e ± e + e~ are unsuppressed by any inverse powers of m R — 
they are suppressed only by small Yukawa couplings — since they may be mediated by 
neutral scalar particles. 

We now want to enforce maximal atmospheric neutrino mixing through an M. u like in 
Eq. (J3*|). We do this by imposing the following generalized CP transformation [7| I12j: 

v Ra -> iS al 3j CP Rl3 , a R -> iS a /3j°Cf3 R , (33) 

01,2 ~> 01,2 > 03 -> -03 ' 

This CP symmetry makes f e real and f^ = f*, while (72^ = 9\r an d #3^ = ~ 9 3r - Without 
loss of generality we assume that v\ is real and positive. Then we find 

M* D = SM D S and M R = SM R S , (34) 

where the second relation follows from the CP invariance of Cm- Therefore too, 
fulfills Eq. ()iup . just as we wanted. 

Let us now consider the masses of the \x and r leptons. Those masses are given by 

"V = \92^2 + 93^3 1 , m T = -j= \g 2T v* 2 + g 3r v 3 *\ = \g* 2fJi v* 2 - g^v* 3 \ . (35) 
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With the notation t> 2j 3 
1 



|^2 3| e 1 ^ 2 ' 3 , we obtain 



m. 



#2 M M +93,x \v3\e 



i(l?2-1?3) 



V2 

Therefore, m T 7^ m M requires 



v/2 



i(l? 3 -1? 2 ) 



,i(*2-*3) ^ _gi(*3-^2) 



Now let us check the case of CP conservation. There we have ^3] 

^2=^2, ^3 = ~ V 3 > 

and therefore 



CP conservation 



; <(tfa-*s) _ _y_ 



(36) 

(37) 

(38) 
(39) 



Thus, if CP is conserved we have = m T ; CP violation is necessary for m M 7^ m T (for 
an earlier model of this type see Ref. [Ti]). 

We have thus shown that, provided CP is spontaneously broken, we are able to obtain 
rrifj, 7^ m T while M. v satisfies Eq. (110)1 . We thus have a model with maximal atmospheric 
mixing but a free \U e ^\. 



4 Obtaining <^ m T 

Let us again consider Eq. ()35j1 . Since m M and m T are both given by essentially the same 
VEVs and Yukawa couplings, it seems natural to expect that they will be of the same 
order of magnitude, even if spontaneous CP breaking ensures that they are different. 
However, in reality one has <C m T . This strong inequality requires in our model the 
almost complete cancellation of two different products of a VEV and a Yukawa coupling. 

In general, if the muon mass is generated by the Yukawa coupling 0t UrD^, where 
is some Higgs doublet, then there is a natural explanation for the smallness of m M : one 
just introduces the symmetry \ir — > —fiR, — > —<f). This symmetry in general restricts 
the Higgs potential in such a way that it allows a vacuum with (0 \<p°\ 0) = 0. The muon 
mass then turns out to be zero. If one now lets the symmetry — > — be softly broken by 
some terms of dimension two in the Higgs potential (0V> where 0' is some other Higgs 
doublet), then we obtain a technically natural explanation for the smallness of (0 \<p°\ 0) 
and thus of m M . Of course, this only works if does not have any Yukawa couplings 
besides the one to the /x, which is not the case in the Standard Model since that model 
only contains one Higgs doublet. 

This idea may be implemented in the context of our model in the previous section. 
Let us introduce the extra symmetry |5, 

K : fi R -> -fi R , 02 ^ 03 (40) 

into that model, then one obtains = —93^ and 

™>ia = \ v 2 - v 3 \ , m T = — j=r \v 2 + v 3 \ . (41) 
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On the other hand, the symmetry K will also restrict the scalar potential, and this in such 
a way that there will in general be a range of parameters of the potential for which the 
vacuum leaves that symmetry unbroken, i.e. v 2 — v^. This immediately leads to m M = 0. 
(Notice that v 2 = ^3 constitutes a maximal spontaneous breaking of our non-standard 
CP symmetry, cf. Eq. ()39j) .) In order to obtain a non-zero but small m M it is now enough 
to introduce into the potential terms which break K softly. In our model, considering its 
other symmetries Tj 2 and CP, there are two such terms: 

l*el (0202 - 0303) + «A*s2 (0203 - 0^) ■ (42) 

The constants /x s i and fi s2 are real. A detailed analysis of the potential is outside the 
scope of this paper [THj, but it is intuitive to expect that, provided /x s i and fi s2 are small, 
\v 2 — will also be kept small and thus m M <C m T will prevail. The crucial point is 
that this inequality is now protected by a softly-broken symmetry, and it is therefore 
technically natural. 



5 Leptogenesis 

Seesaw models offer an attractive possibility for explaining the observed baryon asymme- 
try of the universe by leptogenesis ^Hl EZ] • In order to analyse leptogenesis one works 
in the weak basis where the Majorana mass matrix Mr of the right-handed neutrinos 
is diagonal, with matrix elements Mi, M 2 , and M 3 . The Dirac mass matrix Mb is not 
diagonal in that weak basis. Consider the Hermitian matrix R = MnM.\ 
basis. Then, the CP- violating asymmetry relevant for leptogenesis is 



MdM d in that weak 



3Mi 



IQnv 2 R 



ii 



3 

E 

i=2 



Im 



(Ruf 



Mi 



(43) 



where for simplicity we have assumed Mi -C M 2 , M3 . 

In our model, in the weak basis where Md is diagonal, which is the one that we have 
used before, one has M D = M D = diag(/ e , f^, /*) with f e real, while M^ = SM R S. 
Now, we know from Section 2 that a matrix Mr satisfying M^ = SMrS is diagonalized 
by a unitary matrix V of the Harrison-Scott type, i.e. 



V T M R V = M r = diag {Mi , M 2 , M 3 ) with V 



Pi P2 P3 
Ql <?2 <?3 
Ql (1*2 Q*3 



(44) 



the pj being real. We now move to the weak basis where Mr is diagonal (Mr = Mr) and 
Md is not. In that weak basis Md = V T Mp. Thus, 



/ Pi Qi Q* 

R = V T M D M D V* = p 2 q 2 q* 2 ] diag (/ e 2 , \f, 

V P3 qs ql 



2 

Ml -I 



If, 



( Pi P2 P3 
Q*l <1*2 Ql 

V Qi 12 qs 



(45) 



It is clear that the matrix R is not only Hermitian but, as a matter of fact, real. Therefore, 
in our model leptogenesis is not possible. The root of this fact can be traced directly to 



9 



the existence of a CP symmetry in our model, which is spontaneously broken only at 
the weak scale, i.e. much below the scale at which a net lepton number is supposed to 
be generated. Since CP is unbroken at that super-high scale, e is necessarily zero and 
leptogenesis cannot proceed. 

6 Conclusions 

In this paper we have discussed the mass matrix of Eq. Q, originally found by Babu, 
Ma, and Valle [Hj in the context of a model based on the group A4 and on softly broken 
supersymmetry with additional heavy charged-lepton singlets, an enlarged scalar sector, 
and the seesaw mechanism. In that model, the relations a = b and r = s = hold 
at the seesaw scale and the full mass matrix of Eq. Q arises at the weak scale after 
the renormalization-group evolution of M. v . (Note, however, that subsequently a much 
simpler non-supersymmetric A4 model was proposed by E. Ma [Hj, where the radiative 
corrections to a = b and r = s = are generated by an A<± triplet of charged scalars.) 

Firstly, we have shown that the mass matrix (jHJ) can be characterized by the algebraic 
relation in Eq. (|10|). If we consider all the parameters of the matrix Q as independent, 
it follows readily from this characterization that atmospheric neutrino mixing is maximal 
and that either #13 = (and then the CP phase in lepton mixing is physically meaning- 
less), or #13 is arbitrary and — in the phase convention of Eq. ()20j) — the CP phase is given 
by 7r/2 — see Eq. (|2()jl : this is the more general case. Moreover, the neutrino mass matrix 
of Eq. (JSJ) fixes neither the neutrino masses nor the solar mixing angle. 1 

Secondly, we have derived this mass matrix in the context of a model based on the 
lepton sector of the Standard Model with three right-handed neutrinos, the seesaw mech- 
anism, and three — instead of one — Higgs doublets. We have constructed our model in 
two steps: 

1. Inspired by Refs. j3J|H], we have imposed the three U(l)L a symmetries associated 
with the family lepton numbers, which are softly broken by the Cm of Eq. (|31|). 

2. As suggested by the relation ()10|) and by Ref. j7j, we have imposed the non-standard 
CP symmetry of Eq. (J33|) in order to get the neutrino mass matrix of Eq. (J3J). 

In this way, we have obtained a renormalizable model where lepton mixing arises solely 
from the Majorana mass matrix Mr of the heavy neutrino singlets and where 7^ m T 
is a consequence of the spontaneous breaking of the non-standard CP symmetry. 

Our model contains a symmetry which is spontaneously broken at the Fermi scale 
(see the discussion of the Yukawa Lagrangian (j3~Uj) ). and this may lead to a cosmo logical 
problem through the formation of domain walls at that scale. An additional symmetry 
(see Eq. ()40j0 solves the problem of <C m T in a technically natural way. 

We stress that in our model the neutrino masses are completely free, contrary to what 
happens in the A 4 models which predict neutrinos to be approximately degenerate. We 
also stress that in our model the M. v of Eq. (jSJ) holds at the seesaw scale and its form 

1 As hinted at above, in the A4 models of Ref. [H] the parameters of M. v are not completely independent 
and, therefore, statements about the neutrino mass spectrum and #13 can be made. 
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will be slightly changed by the renormalizat ion-group (RG) evolution down to the Fermi 
scale, while in the A 4 models Eq. (J3J) holds precisely after the RG evolution. Our realiza- 
tion of the mass matrix (jSJ) is an interesting illustration of the fact that exact maximal 
atmospheric neutrino mixing and a non-zero mixing angle #13 can coexist, enforced by a 
symmetry. This was not the case in the models of Refs. ^ E] , where the mass matrix (0) 
was obtained. 

Acknowledgement The work of L.L. was supported by the Portuguese Fundagao para 
a Ciencia e a Tecnologia under the contract CFIF-Plurianual. 
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